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Abstract 

We construct a two-parameter family of Backlund transformations for the trigono- 
metric classical Gaudin magnet. The approach follows closely the one introduced 
by E.Sklyanin and V.Kuznetsov (1998,1999) in a number of seminal papers, and 
takes advantage of the intimate relation between the trigonometric and the ra- 
tional case. As in the paper by A. Hone, V.Kuznetsov and one of the authors 
(O.R.) (2001) the Backlund transformations are presented as explicit symplectic 
maps, starting from their Lax representation. The (expected) connection with 
the xxz Heisenberg chain is established and the rational {xxx) case is recovered in 
a suitable limit. It is shown how to obtain a "physical" transformation mapping 
real variables into real variables. The interpolating Hamiltonian flow is derived 
and some numerical iterations of the map are presented. 



KEYWORDS: Backlund Transformations, Integrable maps, Gaudin systems. Lax rep- 
resentation, r-matrix. 
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1 Introduction 



Backlund transformations are nowadays a widespread useful tool related to the theory 
of nonlinear differential equations. The first historical evidence of their mathematical 
significance was given by Bianchi |3] and Backlund [2] on their works on surfaces of con- 
stant curvature. A simple approach to understand their importance can be to regard 
them as a mechanism allowing to endow a given nonlinear differential equation with a 
nonlinear superposition principle yielding a set of solutions through a merely algebraic 
procedure [19] , [1] , [12] . Backlund transformations are indeed parametric families of dif- 
ference equations encoding the whole set of symmetries of a given integrable dynamical 
system. For finite-dimensional integrable systems the technique of Backlund transfor- 
mations leads to the construction of integrable Poisson maps that discretize a family of 
continuous flows [27] , [25] , [2l] , [22] , [10] , [9] . Actually in the last two decades numerous 
results have appeared in the field of exact discretization of many-body integrable sys- 
tems employing the Backlund transformations tools [IZ] , [21] , [16] , [9] , [10] , [15] , [22] . For 
the rational Gaudin model such discretization has been obtained ten years ago in [S]; 
afterwards, these results have been used for constructing an integrable discretization of 
classical dynamical systems (as the Lagrange top) connected to Gaudin model through 
Inonu-Wigner contractions [T3].|ll].|14]. 

The aim of the present work is to construct Backlund transformations for the Gaudin 
model in the partially anisotropic (xxz) case, i.e. for the trigonometric Gaudin model. 
We point out that partial results on this issue have already been given in [T8] . 

The paper is organized as follows. 

In Section we review the main features of the trigonometric Gaudin model from 
the point of view of its integrability structure. For the sake of completeness, in Section 
([3]) we briefly recall the preliminary results on Backlund Transformations (BTs) for 
trigonometric Gaudin given in [18] . In Section (|4]) the explicit form of BTs is given; it 
is shown that they are indeed a trigonometric generalization of the rational ones (see 
[S]) which can be recovered in a suitable ("small angle" ) limit. The simplecticity of 
the transformations is also discussed in the same Section and the proof allows us to 
elucidate the (expected) link between the Darboux-dressing matrix and the elementary 
Lax matrix for the xxz Heisenberg magnet on the lattice. We end the Section by 
mentioning an open question, namely the construction of an explicit generating function 
for these Backlund transformations. In Section we will show how our map can lead, 
with an appropriate choice of Backlund parameters, to physical transformations, i.e. 
transformations from real variables to real variables. In the last Section we show how a 
suitable continuous limit yields the interpolating Hamiltonian flow and finally present 
numerical examples of iteration of the map. 

2 Gaudin magnet in the trigonometric case 

For a full account of the integrability structure of the classical and quantum Gaudin 
model we refer the reader to the fundamental contributions by Semenov-Tian-Shanski 
[26] and Babelon-Bernard- Talon [1]. In this section we briefly recall the main features 



2 



of the trigonometric Gaudin magnet. 

The Lax matrix of the model is given by the expression: 



L(A) 



M(A) 5(A) 
I C(A) -Ai\) 



N N - 



TV 



(2) 

In (pP) and ([2]) A G C is the spectral parameter, A^ are arbitrary real parameters of 
the model, while (s^, sj, s^) , j = 1, . . . , A^, are the dynamical variables of the system 
obeying to ®^ sl{2) algebra, i.e. 



-2i5jksl, 



(3) 



By fixing the N Casimirs (s^)^ + s'^ s'J = s'j one obtains a symplectic manifold given 
by the direct sum of the correspondent two-spheres. 

Reformulating the Poisson structure in terms of the r-matrix formalism amounts to 
state that the Lax matrix satisfies the linear r-matrix Poisson algebra (see again 



{L(A) 11,11® L(/i)} = [n(A-/i),L(A)®ll+ll®L(/i)]^ 
where rt{X) stands for the trigonometric r matrix [3]: 



(4) 



sin(A) 



/ cos(A) \ 

10 

10 

\ cos(A) / 



(5) 



Equation (|4]) entails the following Poisson brackets for the functions ([2]): 

{A(A), A(/i)} = {5(A), i?(/i)} = {C(A), C(/i)} = 0, 

.cos(A -/i)5(/i) -5(A) 



{A(A),i?(/i)} = z 
{A{X),C{fi)} = z 
{5(A), C(^)} 



sin(A — /i) 
C(A) - cos(A - /i)C(/i) 
sin(A — yu) 
^.2(A(/i)-A(A)) 



(6) 



sin(A — /i) 

The determinant of the Lax matrix is the generating function of the integrals of motion: 



N . 

- det(L) = A\X) + 5(A)C(A) = ^ I 

i=i ^ 



sin (A — Aj 



Hi cot(A — Aj) 



(7) 
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where the Hamiltonians if, are of the form: 



N 



2 COs(Ai - Xk)slsl + S+Sfc +3^3 



k 



Note that only A^— 1 among these Hamiltonians are independent, because of Hi = 0. 
Another integral is given by Hq, the projection of the total spin on the z axis: 

^o = E4 = ^' (9) 

The Hamiltonians Hi are in involution for the Poisson bracket (j3]): 

{H,,Hj} = t,j = 0,...,N-l (10) 
The corresponding Hamiltonian flows are then given by: 

In the XXX model a remarkable Hamiltonian is found by taking a linear combination 
of the integrals corresponding to ([8]) in the rational case [6]. It describes a mean field 
spin-spin interaction: 

Where the notation for the bold symbol Sj is Sj = {3}, s^, 3^) with 3^ = 3} + is^ and 
3~ = 3] — 131. The natural trigonometric generalization of this Hamiltonian can be 
found by taking the linear combination of [HI 

sin(2Aj) 
2 ' 

4 = 1 

giving 

1 ^ 

•Ht = - cos(A, + A,) {3]3] + s^s^' + cos(Ai - Aj)sf s^^) (12) 
*7^i 



TV 

2 ^* ■ 



3 A first approach to Darboux-dressing matrix 

In this Section, for the sake of completeness, we recall the results already appeared in 
[T8] . The leading observation is that by performing the "uniformization" mapping: 

the TV-sites trigonometric Lax matrix takes a rational form in z that corresponds to 
the 2iV-sites rational Lax matrix plus an additional reflection symmetry (see also [7j); 
in fact, by performing the substitution ([3]), the Lax matrix ([T]) becomes: 
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Liz) =^J' + f2 - as^as ] , (13) 



where is the Pauh matrix diag{l, —1) and the matrices L{, j = 1, . . . , N , are given 
by: 



L{ = iz 



So, equation f lT3|) entails the following involution on L{z): 

Liz) = a3L{-z)a3 (14) 

Constructing a Backlund transformation for the Trigonometric Gaudin System (TGS) 
amounts to build up a Poisson map for the field variables of the model ([2]) such that the 
integrals of motion (|8]) are preserved. At the level of Lax matrices, this transformation 
is usually seeked as a similarity transformation between an old, or "undressed". Lax 
matrix L, and a new, or "dressed" one, say L: 

L{z) D{z)L{z)D-\z) = L{z) (15) 

But L and L have to enjoy the same reflection symmetry f ll4p too: to preserve this 
involution the Darboux dressing matrix D has to share with L the property f lT^ : 
the elementary dressing matrix D is then obtained by requiring the existence of only 
one pair of opposite poles for D in the complex plane of the spectral parameter. We 
will show in the next Section that, thanks to this constraint, one recovers the form 
of the Lax matrix for the elementary xxz Heisenberg spin chain: on the other hand, 
this is quite natural if one recalls that for the rational Gaudin model the elementary 
Darboux-dressing matrix is given by the Lax matrix for the elementary xxx Heisenberg 
spin chain [8], [10]. The previous observations lead to the following Darboux matrix: 

D{Z) = Doc + - ^3-^^3 (16) 

z-^ z+^ 

By taking the limit 2; — )■ oo in ( !T6|) it is readily seen that Dao has to be a diagonal 
matrix. In order to ensure that L and L have the same rational structure in z, we 
rewrite equation ( !T5|l in the form: 

L{z)D{z) = D{z)L{z) (17) 

Now it is clear that both sides have the same residues at the poles z = zj, z = (it 
is unnecessary to look at the poles in 2; = —Zj and z = — because of the symmetry 
f|T^ . so that the following set of equations have to be satisfied: 

L?D{z,) = D{z,)L?, (18) 
LiODi = D,L{0- (19) 
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In principle, equations (fT8|) . ( IT9|) yield a Darboux matrix depending both on the old 
(untilded) variables and the new (tilded) ones, implying in turn an implicit relationship 
between the same variables. To get an explicit relationship one has to resort to the so- 
called spectrality property [TU] [H]. To this aim we need to force the determinant of the 
Darboux matrix D{z) to have, besides the pair of poles a.t z = a pair of opposite 
nondynamical zeroes, say ai z = ±77, and to allow the matrix Di to be proportional to 
a projector [18]. Again by symmetry it suffices to consider just one of these zeroes. If 
7] is a zero of detD{z), then -D(?7) is a rank one matrix, possessing a one dimensional 
kernel \K{ri)); the equation (fT7|) : 

L(r/)D(r/) = I)(r/)L(r/) (20) 

entails 

D{v)L{r^)\Kiv)) = 0. (21) 
This equation in turn allows to infer that \K{r])) is an eigenvector for the Lax matrix 

L{rJ)\Kil^))=^^iv)\K{l^)), (22) 

This relations gives a direct link between the parameters appearing in the dressing 
matrix D and the old dynamical variables in L. Because of f[T^ we have another one 
dimensional kernel \K{^)) of Di, such that: 

m\m)) = m\m))- m 

In [18] we have shown how the two spectrality conditions (122|) . fl23l) enable to write D 
in terms of the old dynamical variables and of the two Backlund parameters and rj. 
The explicit expression of the Darboux dressing matrix is given by: 



b ^ bz b . 

bp{mv){e-V^) fc(p(rj)rj-p(g)0 bz{p{Ov-p{ri)0 

z ^ r,i 

In this expression /3 is a global multiplicative factor, inessential with respect to the 
form of the BT, b is an undeterminate parameter that in Section (jl]) we will fix in 
order to recover the form of the Lax matrix for the discrete xxz Heisenberg spin chain. 
The functions p{ri) and p(^) characterize completely the kernels of -D(?7) and -D(0- iii 
fact we have the following formulas [T8] : 



mo) - ( ) mr,)) ^ ( ^j^, ) (25) 

As \K{^)) and \K{ri)) are respectively eigenvectors of L(^) and L{ri), p{^) and p{ri) 
must satisfy: 

Mr,>^^^ (26) 

with A{z), B{z), C{z) given by 1^ and iJ,\z) = A^{z) + B{z)C{z). 
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4 Explicit map and an equivalent approach to Darboux- 
dressing matrix 



The matrix ( l24l) contains just one set of dynamical variables so that the relation ( |T5i) 
gives now an explicit map between the variables (s^, sj, s|) and (s^, sj , s|) . The map 
is easily found by ( |T8l) : it reads: 

~3 _ PiQpiv) (^^ - V^) ii^k^ - V^) PiO^ - W - i'^) Piv)v) SkZk , 



{e - v') iizk' - e) PiOv - pm i^k' - v')) stz. 



A, 

PiOpiv) ae + zk') iv' + zk') - iv' + e) - sv^e^k') (^ra) 



A, 

ivae-z,') iv'-zk') (pio'+piv?)) 



+ 



2bM0piv) ie - v') iizk' - e) PiOv - pm izk' - v')) sizk 



(27b) 



W slzl^iq ^ {{zk^ -V^)p{0^-{zk^ -^^)p{v)vf Sf^^r] ^ 

(27c) 



2 ie - v') {{zk' - v') Pirn - izk' - e) p^) sUk^v 



b^A, 

where is proportional to the determinant of D{zk), i.e. 

Ak = {zl -e)izl-v') (piOv - pm) iPiv)v - PiOO (28) 

Formulas f l27ap . (]27bp . fl27cp define a two-parameter Backlund transformation, the 
parameters being ^ and rj: as we will show in the next section, it is a crucial point to 
have a iwo-parameter family of transformations when looking for a physical map from 
real variables to real variables. As mentioned in the previous Section, we now show 
that indeed, by posing: 

b = i^f^ (29) 

the expression (!24l) of the dressing matrix goes into the expression of the elementary 
Lax matrix for the classical, partially anisotropic, Heisenberg spin chain on the lattice 
0. 
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Obviously two matrices differing only for a global multiplicative factor give rise to the 

1^ 



same similarity transformation. So we omit the term -^rzp in ( 124|) . and, taking into 



account fl29|l . we write for the diagonal part Dd of (|2^ : 

Dd = l{ipiO - Piv)){v - w)l+ip{0+p{v))iv + ^)^3) (30) 
where v{C,,ri) and w{C,,f]) are given by: 



-(e, r/) = ^ - ^ r/) = ^ - ^ = -viv, (31) 



We substitute: 

^ ^ e^^i r/ ^ e^^2 z e'^ (32) 

and take a suitable redefinition of the Backlund parameters to clarify the structure of 
the D matrix: 

^0 = — ^ — = — ^ — (33) 

With these positions it is simple to find that v — w = Aie^^° sin(A — Aq) cos(/i) and 
V + w = 4ze*'^° cos(A — Aq) sin(/i). So, considering equation ( l30l) jointly with the off- 
diagonal part of fl2^ . the dressing matrix can be written as: 



D{\) =a sin(A-Ao)U 
2 sin(/x) 

p(C2)-p(Ci) 



+ 



P(Cl) + P(C2) 

p(Ci)-p(C2) 

1 

-p(Ci)p(C2) 



tan(/i) cos(A — Ao)o"3+ 



(34) 



where a is the global factor 2e*^''(p((^2) ~p(Ci))- Observe that in formula f lM|) . with 
some abuse of notation, p(Ci) (p(C2)) stands of course for p(OI^=e<i (^P('7)lr;=e<2 
A last change of variables allows to identify the dressing matrix with the elementary 
Lax matrix of the classical xxz Heisenberg spin chain on the lattice, and furthermore 
to recover the form of the Darboux matrix for the rational Gaudin model [S] [2D] in the 
limit of small angles. Namely, we introduce two new functions, P and Q, by letting 



p(Ci) = -Q P{C2 



2 sin(/i) 
P 



Q. 



(35) 



Then equation flMj) becomes: 



D{X) = a 



sin(A — Ao — /i) + PQ cos(A — Aq 



P cos(/i) 
sin(A — Ao + /i) — PQ cos(A — Aq 



(36) 

Obviously now we can repeat the argument made before about spectrality; indeed now 
D L , and D L , are rank one matrices. So if il+ and are respectively the 
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kernels of -D(Ao + /i) and -D(Ao — fi) one has again that f2+ and f2_ are eigenvectors of 
L(Ao + /i) and L{\q — jj) with eigenvalues 7+ and 7_ where 



A=Ao±/^ 



7± = 7(A) 
and we have set ([7]) 

72(A) = A2(A) + B{\)C{\) = -det(L(A)) 
The two kernels are given by: 



-Q 



2 sin(^) - PQ 



(37) 



(3^ 



and the eigenvectors relations yields the following expression of P and Q in terms of 
the old variables only: 



Q = Q(Ao + /i) 



A(A)-7(A) 



i?(A) 



1 g(Ao + /i) -g(Ao-/i) 



A=Ao+A» -P 



2 sin(/i) 



(39) 



The explicit map can be found by equating the residues at the poles A = A^ in ( 1T7|) . 
that is by the relation: 

hDk = DuLk (40) 

where 



st -si " 



Dk = D{\ = Afc) 



(41) 



or by performing the needed changes of variables in fl27ap . (]27bp . ( 12 7c p . Anyway now 
the map reads: 



2 cos\fi) - (cos2(/i) + cos2((5g))(l - 2PQ sin(/i) + P^Q"' 

, Pcos(/i)(sin(,5^)-Pgcos(^g)) ^ 

+ A. + 

Q cos(/i)(2 sin(/i) - Pg)(sin(5^) + PQ cos(5^)) 



(42a) 



sm{6l) - PQ cos{5^)f ^ (g2 cos2(/i)(2 sin(/i) - PQ)f 



A, 



A, 



2Q cos(;u)(2 sin(^) - Pg) (sin(5^) - PQ cos{5^)) 



A, 



_ _ (sin((5^) + Pg cos(^g))2 _ _ P2 cos2(/i) + 
A, A, + 

2P cos(;u)(sin(5^) + PQ cos(5^)) 3 



(42b) 



(42c) 
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where for typesetting brevity we have put: 

^0 = Afc - Ao / .o\ 

6i = Xk-\o±fi ^ ' 

and we have denoted by the determinant of D{Xk), that is: 

Afc := sin(Afc - Ao - yu) sin(Afc - Ao + - 2PQ sin(/i) + P^Q^) 

At this point we can show that for "small" Aq and fi one obtains, at first order, the 
Backlund for the rational Gaudin model, independently found by Sklyanin [20] on one 
hand and Hone, Kuznetsov and Ragnisco [8] on the other, as the composition of two 
one-parameter Backlunds. So let us take Aq — )■ hXo, ^ hfj, and X hX where h is 
the expansion parameter. One has: 

cot(A - Afc) = \ . + 0{h) ■ , . = . + 0{h), 

h{X - Xk) sm(A - Afc) h{X - Xk) 

so that Q = + 0{h'^), where the superscript r stands for "rational" . Thus, 
coincides with the variable q that one finds in the rational case [8] . For the variable P 
one has: 

P = h{p' + 0{h^)) where / 



g''(Ao + /u) - g''(Ao - /x) 

Taking into account these expressions, it is straightforward to see that the matrix (l36l) 
has the expansion: 

D{X) = hD'{X) + 0{h^) (44) 

where 

^^(A)=f' -4°"^^?;^'^ , , / V (45) 

^ ' \ q' [2fi — p^q^) A — Ao + /i— / 

The limit of "small angles" in f l27ap . ( I27bl) . fl27cp obviously leads to the rational map 

of El. 



4.1 Symplecticity 

In this subsection we face the question of the simplecticity of our map; the corre- 
spondence with the rational Backlund in the limit of "small angles" shows that the 
transformations are surely canonical in this limit. Indeed, as our map is explicit, we 
could check by brute-foce calculations whether the Poisson structure ([3]) is preserved 
by tilded variables. However we will follow a finer argument due to Sklyanin [2T] . 
Suppose that -D(A) obeys the quadratic Poisson bracket, that is 

{D\X), D\r)} = [n{X - r), D\X) ® D\r)] (46) 

where as usually = D ® )[, = 11 ® -D. Consider the relation 

L(A)^(A - Ao) = D{X - Ao)L(A) (47) 
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in an extended phase space, where the entries of D Poisson commutes with those of L. 
Note that in (H71) we have used tilded variables also for -D(A) (in its l.h.s.) because (l47j) 
is indeed the Backlund transformation in this extended phase space, whose coordinates 
are (sj,s^,P, Q), so that we have also a P and a Q. The key observation is that if 
both L and D have the same Poisson structure, given by equation (H^ . then this 
property holds true for LD and DL as well, because in this extended space the entries 
of D Poisson commute with the entries of L. This means that the transformation fHTj) 
defines a "canonical" transformation. Sklyanin showed [21] that if one now restricts the 
variables on the constraint manifold P = P and Q = Q the symplecticity is preserved; 
however this constraint leads to a dependence of P and Q on the entries of L, that for 
consistency must be the same as the one given by the equation ( H7I) on this constrained 
manifold. But there (147|) is just given by the usual BT: 

L(A)D(A-Ao) = D(A-Ao)L(A) 

so that the map preserves the spectrum of -^^(A) and is canonical. What remains to 
show is that indeed (146|) is fullfilled by our -D(A). Obviously -D(A) cannot have this 
Poisson structure for any Poisson bracket between P and Q. In the rational case the 
Darboux matrix has the Poisson structure imposed by the rational r-matrix provided 
P and Q are canonically conjugated in the extended space [21] (and this is why they 
were called P and Q); in the trigonometric case P and Q are no longer canonically 
conjugated but obviously one recovers this property at order h in the "small angle" 
limit. 

First note that -D(A) can be conveniently written as: 

D{\) = a cos(/x) sin(A)ll + a cos(A)(T3 + ^ q ^ (48) 
where the coefficients a, b, c are given by: 

Pg- sm(/i) ^ 6 = P, c = 2Qsm{^^)-PQ' (49) 
cos(/x) 

Inserting in (HSl) we have the following constraints: 

{a,aa} = =^ a = a{PQ) (50) 

{a,ab} = -a'ab =^ p} = aP^l^M^-^ (51) 

cos(/i) 

2 r ^^ ^Sin(/i) — PQ 

{a,ac\ = aac =^ {a,Q \ = —aQ — — (52) 

cos(/i) 

All remaining relations, namely 

{ab, ac} = 2c?a {aa, (xh\ = a^b {aa, ac\ = —o?c (53) 

give the same constraint, i.e.: 

Xn PI l + P^g^-2Pgsin(/x) 
cos(/i) 
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This expression can be used to find, after a simple integration, 



a{PQ) = , 

^(l + p2g2_2Pgsin(/i)) 

so tliat tlie Darboux matrix (136|) is fixed (up to the constant multiphcative factor k). As 
previously pointed out, it turns out that the Darboux-dressing matrix (136|) is formally 
equivalent to the elementary Lax matrix for the classical xxz Heisenberg spin chain 
on the lattice [5]. Moreover it has also the same (quadratic) Poisson bracket. This 
suggests that indeed -D(A) can be recast in the form (see [5]): 

D{X) = ^ol + -^{yiai + ^2^2 + cos(A)Jl3a3) (55) 
sm(^A j 

where the cxj are the Pauli matrices and the variables satisfies the following Poisson 
bracket (0): 

where {i,j,k) is a cyclic permutation of (1,2,3) and Jjk is antisymmetric with Ju = 
0, Ji3 = J23 = 1. Indeed it is straightforward to show that the link between the two 
representations ( 148|) and ( l55l) . up to the factor cos(/i) sin(A) that does not affect neither 
( ITTl) nor the Poisson bracket ( l46l) . is given by : 

a = ^o -y(& + c) = ^i -(fe-c) = ^2 -ma = ^3 (57) 



and the Poisson brackets ( l50l) . ( 1511) . (1521) . ( 1531) correspond to those given in ( 1561) . 
An open question regards the generating function of our BT. So far we have not been 
able to write it in any closed form; in our opinion the question is harder than in the 
rational case (where the generating function is known from |8]): in fact the rational 
map corresponding to (I27ap . ( ]27b[) . (I27cp can be written as the composition of two 
simpler one-parameter Backlund transformations, and this entails the same property 
to hold for the generating function; in the trigonometric case a factorization of the 
Backlund transformations cannot preserve the symmetry (fT4|) : so probably one should 
look for symmetry-violating generating functions such that their composition enables 
symmetry to be restored. 



5 Physical Backlund transformations 

The transformations we have found do not map, in general, real variables into real 
variables. A sufficient condition to ensure this property is given by: 



Ci = C2 (58) 



which amounts te require that Aq and /x in (I42ap . ( ]42bl) . (I42cp be, respectively, real and 
imaginary numbers. 
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Indeed we claim that, if ( 158|) holds, starting from a physical solution of the dynamical 
equations, we can find a new physical solution with two real parameters. Let us prove 
the assertion. Backlund transformation are obtained by (HOl) : starting from a real 
solution means starting from an Hermitian L^. Thus, if the transformed matrix has 
to be Hermitian too, the Darboux matrix has to be proportional to an unitary matrix. 
We will show that this is the case by choosing d = (2 and 7(Ci) = — 7(C2) (7 is the 
function defined in fl37j) ). Note that the condition on the 7's specifies their relative sign 
(the sheet on the Riemann surface), inessential for the spectrality property. Hereafter 
we assume the parameter fi, defined in ( l32l) . to be purely imaginary = ie, so that: 

Ci = Ao + le (Ao, e) G (59) 

The Darboux matrix at A = A^ can be rewritten as: 

^ / sm{vk - ie) + PQ cos{vk) P cosh(e) \ , . 

~ V Q cosh(e) (2z sinh(e) - PQ) sm{vk + le) - PQ cos{vk) J ^ ' 

where = Xk — Aq (we are assuming that the parameters A^ of the model are real). 
We recall that in (lUOj) : 

^_^..^_ -4(Cl)-7(Cl) _ C(Ci) . o _ 2z sinh(6) 

^ ^^^'^ BiCi) A(Ci) + 7(Ci)' g(Ci)-Q(Ci)' ^ ^ 

Furthermore it is a simple matter to show that 

A(Ci) = A(Ci); i?(Ci) = C'(Ci); C(Ci) = fi(Ci)- (62) 

If the off-diagonal terms of D^dI has to be zero, then the following equation has to be 
fuimiled: 

P(sin(t;fc - le) - PQ cos{vk)) = Q{2i sinh(e) + PQ)ism{vk - le) + PQ cos{vk)) (63) 
Using relations fl6T]) and rearranging the terms, the previous equation becomes: 



11, 11 

' cosh(e) sin(wfc) + z( - + - - ) cos(t'fc) sinh(e) 



ViCi) QiCi)' ' ' ^Q(Ci) Q(Ci)' ' ' (64) 

= (Q(Ci) - Q(Ci)) cosh(e) sm{vk) + i cos{vk) sinh(e)(Q(Ci) - g(Ci)) 



Note that the relations fl62l) gives 7^(Ci) = 7^(Ci)) implying that 7^ (A) is a real function 
of its complex argument, consistently with the expansion ([7]). 
The choice: 

7(Ci) = -7(Ci) (65) 

entails: 

With this constraint the equation (164|1 holds too. Moreover (165|) makes the diagonal 
terms in DkD\ equal. This shows that, under the given assumptions, Dk is an unitary 
matrix. 
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6 Interpolating Hamiltonian flow 

The Backlund transformation can be seen as a time discretization of a one-parameter 
(Ao) family of hamiltonian flows with the difference — (i) = 2e playing the role of 
the time-step. To clarify this point, let us take the limit e — 0. 
We have: 

Q = ^^^^ + Oie)^Qo + Oie) (67) 

P = -te^^ + 0(e2) = lePo + 0(6^) (68) 
7(Ao) 

and for the dressing matrix we can write: 
D{\) = k sin(A - Ao)U+ 

,,V.^cos(A-Ao)(Pogo-l) Po Vnr,2^ (69) 

^ go(2 - PoQo) cos(A - Ao)(l - PoQo) J ^ ^ 

Reorganizing the terms with the help of Pq and Qq given in the equations fl57|) and 
we arrive at the expression: 

D{X) = k sin(A - Ao)ll+ 
tek A(Ao) cos(A - Ao) 5(Ao) ^ + 0(6^) ^^^^ 



7(Ao) V C'(Ao) -A(Ao)cos(A-Ao) 

It is now straightforward to show that in the limit e — the equation of the map 
LD = DL turns into the Lax equation for a continuous flow: 

L(A) = [L(A),M(A,Ao)] (71) 
where the time derivative is defined as: 

L = lim (72) 
e^o e 



and the matrix M(A, Aq) has the form 



A(Ao)cot(A-Ao) -Ig^ 

7(Ao)l, -g^ -A(Ao)cot(A-Ao) 



(73) 



The system (1711) can be cast in Hamiltonian form: 

L(A) = {H(Ao),L(A)} (74) 
with the Hamilton's function given by: 



^(Ao) = 7(Ao) = VA2(Ao) + 5(Ao)C(Ao) (75) 



Quite remarkably, but not surprisingly, the Hamiltonian (1751) characterizing the inter- 
polating flow is (the square root of) the generating function of the whole set of 
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conserved quantities. By choosing the parameter Aq to be equal to any of the poles 
(Aj) of the Lax matrix, the map leads to different maps {i?T'-*)}j=i..Ar, where BT^^^ 
discretizes the flow corresponding to the Hamiltonian Hi, given by equation ([8]). Any 
other integrable map for the trigonometric Gaudin model can be, in principle, written 
in terms of the maps {BT^^^}i=i j\i. 

More explicitely, by posing Aq = 5 + Aj and taking the limit 6^0, the Hamilton's 
function flTSjl gives: 

7(Ao) = y + ^ + 0(5) (76) 
2si 

and the equations of motion take the form: 

L(A) = -^{i7„L(A)} (77) 

Accordingly, the interpolating flow encompasses all the commuting flows of the system, 
so that the Backlund transformations turn out to be an exact time- discretizations of 
such interpolating flow. 



6.1 Numerics 




Figure 1: input parameters: Si = 2+z, = 2—i, si = —2, = 50+40z, S2 = 50— 40z, 
si = 70, Ai = 7r/110, A2 = 77r/3, Aq = 0.1, /x = -0.002i 



The figures report an example of iteration of the map fl42al) . fl42bp . fl42cp . For 



simplicity we take N = 2. The computations shows the first 1500 iterations: the 
plotted variables are the physical ones {sf, s\, sf). Only one of the two spins is shown, 
namely that labeled by the subscript "1". The figures are obtained by a "'"'^Maple 
code. 
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Figure 2: input parameters: = 0.2 + lOi, = 0.2 — lOi, sf = — 1, = 10 — 30i, 
s" = 10 + 30i, si = 100, Ai = TT, As = Ttt/S, Aq = 0.1, /i = -0.004i 
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